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Abstract 

We review the Cosmology and Physics underlying Primordial Nucleosynthe- 
sis and survey current observational data in order to compare the predictions 
of Big Bang Nucleosynthesis with the inferred primordial abundances. From 
this comparison we report on the status of the consistency of the standard hot 
big bang model, we constrain the universal density of baryons (nucleons), and 
we set limits to the numbers and/or effective interactions of hypothetical new 
"light" particles (equivalent massless neutrinos). 
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1 Introduction 



At present, the Universe is observed to be expanding jjj and filled with radiation 
U - [| which is very cold today (T = 2.728 K g). If the evolution of such a 
Universe is traced back in time to earlier epochs which were hotter and denser, 
the early Universe is a Primordial Nuclear Reactor during its first 20 minutes 
(psIOOO sec). As the early Universe expands and cools, nuclear reactions are pre- 
maturely aborted before the heavier elements can be synthesized. Only the light 
nuclides deuterium (D), hclium-3 ( 3 He), helium-4 ( 4 He), and lithium-7 ( 7 Li) can 
be synthesized in abundances comparable to those observed (or, observable!) in 
a variety of astrophysical sites (e.g., stars; cool, neutral gas; hot, ionized gas). 
Since the relative abundances of the primordially-produced nuclides depend on 
the density of nucleons (baryons) and on the early-Universe expansion rate, a 
comparison of the predicted and observed abundances provides a key test of 
the standard model of cosmology, as well as an indirect "measurement" of the 
baryon density of the Universe which is equally sensitive to dark and luminous 
baryons (i.e., Is the early-Universe nucleon abundance consistent with that in- 
ferred today from non-BBN data?), and offers a unique probe of hypothetical 
new particles (beyond the standard model) whose presence would have altered 
the expansion rate of the early Universe (hence changing the time available for 
element synthesis) . As one of the pillars of the standard model of Cosmology, 
BBN opens a unique window on the Universe. 

In this review, dedicated to the memory of our friend and colleague Dave 
Schramm, we review the basic physics and cosmology relevant to the calculation 
of the primordial yields, both in the standard model and in simple extensions 
of the standard model. Then we compare the current predictions, based on 
up-to-date nuclear and weak interaction rates, with the primordial abundances 
inferred from observational data obtained from a variety of astrophysical sites 
using a variety of astronomical techniques. Since the BBN-prediction part is 
relatively simple and straightforward, it is the data which lies at the core of 
such comparisons. The good news is that the key nuclides are observed in a 
variety of objects using very different techniques, thus minimizing correlated 
systematic errors in the abundance determinations. Additional good news is 
that modern telescopes and detectors have provided high quality data whose 
statistical errors have been shrinking dramatically. The bad news is that, for 
most abundance determinations, the accuracy is now limited by our ignorance 
of possible systematic errors which are often difficult to quantify using extant 
data alone. Therefore, a large part of this review is devoted to the data and 
our assessment of the uncertainties. In this, we strive to err on the side of 
caution. When mutually contradictory data appears (as it does for primordial 
deuterium) we will explore the consequences of each option, letting the reader 
draw his/her own conclusion. Given the rapid pace of observational cosmology 
at present, the quantitative abundances derived from current data are likely 
ephemeral. However, it is our hope that our discussion here will set the stage 
for any changes new data will provide. 
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2 Primordial Nucleosynthesis 



All that is needed to predict the primordial abundances of the light elements 
within the context of the standard models of cosmology and particle physics 
is the current temperature and expansion rate of the Universe. Then, under 
the assumptions that the Universe is homogeneous and isotropic and that the 
standard model of particle physics is the correct description of the particle 
content of the Universe at temperatures of order a few MeV, the predicted 
primordial abundances of D, 3 He, 4 He, and 7 Li depend only on the baryon 
density. That is, the predictions of standard BBN are uniquely determined by 
one parameter, 77, the baryon-to-photon ratio: 7710 = 273f2B^ 2 (^b is the ratio 
of the baryon density to the critical density and the Hubble parameter is Ho = 
lOOh km/s/Mpc; 7710 = I0 w r]). 

The primordial yields of light elements are determined by a competition 
between the expansion rate of the Universe, the rates of the weak interactions 
that interconvert neutrons and protons, and the rates of the nuclear reactions 
that build up the complex nuclei. Neglecting the contributions of curvature and 
the cosmological constant, which are small in the early Universe, the expansion 
rate is determined by the Friedmann equation: 



where R is the scale factor. For standard BBN the energy density, p, at the 
time nucleosynthesis begins (about 1 second after the Big Bang) is described by 
the standard model of particle physics 



where p 7 , p e , and p v are the energy density of photons, electrons and positrons, 
and massless neutrinos and anti- neutrinos (one species), respectively, and N v is 
the equivalent number of massless neutrino species which, in standard BBN, is 
exactly 3. 

At high temperatures, neutrons and protons can interconvert via weak in- 
teractions: n + e + <-» p + v e , n + v e <-» p + e~, and n <-» p + e~ + v e . As 
long as the interconversion rate of neutrons and protons is faster than the ex- 
pansion rate, the neutron-to-proton ratio tracks its equilibrium value, expo- 
nentially decreasing with temperature. This condition holds for temperatures 
T > 1 MeV as can be seen from a comparison of estimates of the weak rates 
r wk = n{av)~ 0(1CT 2 )T 5 /M^, and the expansion rate H = (8nG N p/3) 1/2 ~ 
5.4 T 2 /Mp, where Mw and Mp are the electroweak and Planck masses, respec- 
tively^]. Once the interconversion rate becomes less than the expansion rate, 
n/p effectively "freezes-out" (at about 1/6), thereafter decreasing slowly due to 
free neutron decay. 

1 A detailed numerical calculation would find that equilibrium is maintained down to 0.8 




(1) 



p = p 1 + p e + N^p, 



(2) 



MeV. 
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Although freeze-out occurs at a temperature below the deuterium binding- 
energy, Eb = 2.2 MeV, the first link in the nucleosynthetic chain, p+n — ► D + 7, 
is not effective since the photo-destruction rate of deuterium (oc n 1 e.~ EB l T ) is 
much larger than the production rate (oc ng) due to the large photon-to baryon 
ratio (> 10 9 ). As soon as deuterium becomes stable against photodissociation 
(~ 80 keV) neutrons are bound up into 4 He with an efficiency of 99.99%, driven 
by the stability of the 4 He nucleus. By this time, n/p has dropped to ~ 1/7, 
and simple counting yields an estimated 4 He mass fraction 

Yp « - 2( " /p) - = 0.25. (3) 
[1 + (n/p)} > 

In addition, the large binding energy of 4 He insures that its primordial abun- 
dance is relatively insensitive to the nuclear reaction rates (which increase with 
increasing baryon density (rf)). D (and 3 He) is burned to get to complex nuclei 
and thus its abundance decreases rapidly with increasing 77, making D the per- 
fect baryometer (see section 3.1). At low 77, 7 Li is destroyed by protons with an 
efficiency that increases with 77, while at relatively high 77, 7 Be (the dominant 
route to 7 Li by subsequent electron capture) is produced more efficiently with 
increasing 77. Hence the "Li valley" in a 7 Li vs. 77 plot. Increasing Coulomb 
barriers and a lack of stable nuclei at A = 5 and 8 cause standard BBN to 
struggle to produce 7 Li and be even less effective beyond that. 

In Figures 1-3, the primordial abundances predicted by standard BBN are 
shown as a function of 77. The width of each curve reflects the 2a uncertainty in 
the predictions that results from a Monte Carlo analysis of uncertainties in the 
neutron lifetime and nuclear reaction rates § . The neutron lifetime was taken 
to be t„ = 887 ±2 seconds^. At 77 = 5 x 10~ 10 the fractional uncertainties due to 
2a experimental errors are 0.4% for Y, 15% for D or 3 He, and 42% for 7 Li if the 
errors in the cross sections from Smith et al. @ are used. We note that recent 
work ||, based on a preliminary reanalysis and update of the relevant reaction 
cross sections, claims smaller uncertainties in D and 7 Li by roughly a factor of 
two. The robustness of the BBN predictions is directly related to the fact that, 
for the most part, the astrophysical S-factors are measured at energies relevant 
to the BBN environment. 

The fractional uncertainty in the predicted mass fraction of 4 He due to ex- 
perimental errors in the reaction rates is, for the range of nucleon density of 
interest, almost entirely due to the uncertainty in the neutron lifetime (which 
translates into an overall uncertainty in the weak interconversion rates). In 
the last few years considerable effort has gone into understanding the theoreti- 
cal uncertainty in the predicted abundance of 4 He due to the treatment of the 
weak interaction rates. The BBN code used for the results presented in Walker 
et al. (WSSOK) g] was basically the 'Wagoner Code' |l§ (updated first by 
Yang with modifications by Walker |I| to include zero and finite- 

temperature radiative corrections and Coulomb corrections to the weak rates 

2 The current world average is r n = 886.7 ± 1.9 Q| leading to predictions indistinguishable 
from those displayed. 
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Figure 1: The predicted 4 He abundance (solid curve) and the 2a theoretical un- 
certainty H . The horizontal lines show the range indicated by the observational 
data. 

(as described in Dicus et al. fill). Subsequent modifications included an up- 
date and enlargement of the nuclear network by Thomas in TSOF [ fl6| . Kernan 
re-examined the codejlT]], updating the TSOF code to include finite nucleon 
mass effects (as described in Seckel |L8| and Lopez, Turner, and Gyuk |Q) and 
finding a relatively large time-step error. He also estimated the uncertainty 
in the 4 He mass fraction due to choice of finite-temperature prescription and 
non-equilibrium effects in the neutrino sector to be ~ 10~ 4 . Kernan's recom- 
mendation was, to a level of accuracy of a few parts in 10 4 , to simply use the 
Y P found in the 'WSSOK Code' and add 0.003. This was adopted in Hata et 
al. H to yield Yp = 0.2467 (theoretical uncertainty of a few parts in 10~ 4 ) at 
77 = 5 x 10~ 10 and r„ = 887 seconds. With several versions of the BBN code 
floating around, no one but the owners (and sometimes not even they) knew how 
the various "Yp -corrections" were handled and no one had built an indepen- 
dent version of the code that contained all these corrections in a self-consistent 
way. Lopez and Turner [po[ have recently done just that. Starting from scratch, 
including all the effects mentioned above, and adding order-a QED corrections 
(as described by Heckler |pH ) and detailed non-equilibrium neutrino effects (as 
described by Dodelson and Turner H they find: Y P = 0.2460 ±0.0002(theory) 
at rj = 5 x 10~ 10 and t„ = 885.4 seconds (if they use r„ = 887 seconds they find 
Yp = 0.2467 (Lopez, Private Communication to G.S.)). Indeed, over the entire 
range 1 < 7710 < 10, the difference in predicted 4 He mass fraction between our 
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Figure 2: The predicted D/H abundance (solid curve) and the 2a theoretical 
uncertainty ||. The horizontal lines show the range indicated by the observa- 
tional data for both the high D/H (upper two lines ) and low D/H (lower two 
lines) . 

code and the Lopez/Turner code is 0.0001 ±0.0001. 

3 From Observations To Primordial Abundances 

To test the standard model it is necessary to confront the predictions of BBN 
with the primordial abundances of the light nuclides which are not "observed" , 
but are inferred from observations. The path from observational data to primor- 
dial abundances is long and twisted and often fraught with peril. In addition 
to the usual statistical and insidious systematic uncertainties, it is necessary to 
forge the connection from "here and now" to "there and then", i.e., to relate 
the derived abundances to their primordial values. It is fortunate that each 
of the key elements is observed in different astrophysical sites using very dif- 
ferent astronomical techniques and that the corrections for chemical evolution 
differ and, even more important, can be minimized. For example, deuterium is 
mainly observed in cool, neutral gas (H i regions) via resonant UV absorption 
from the ground state (Lyman series), while radio telescopes allow helium-3 to 
be studied via the analog of the 21 cm line for 3 He + in regions of hot, ionized 
gas (Hn regions). Helium-4 is probed via emission from its optical recombi- 
nation lines in Hn regions. In contrast, lithium is observed in the absorption 
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Figure 3: The predicted 7 Li abundance (solid curve) and the 2a theoretical un- 
certainty H . The horizontal lines show the range indicated by the observational 
data. 

spectra of hot, low-mass halo stars. With such different sites, with the mix of 
absorption/emission, and with the variety of telescopes involved, the possibil- 
ity of correlated errors biasing the comparison with the predictions of BBN is 
unlikely. This favorable situation extends to the obligatory evolutionary cor- 
rections. For example, although until recently observations of deuterium were 
limited to the solar system and the Galaxy, mandating uncertain corrections 
to infer the pregalactic abundance, the Keck and Hubble Space telescopes have 
begun to open the window to deuterium in high-redshift, low-metallicity, nearly 
primordial regions (Lyman-a clouds). Observations of 4 He in low-metallicity 
(~ 1/50 of solar) extragalactic Hn regions permit the evolutionary correction 
to be reduced to the level of the statistical uncertainties. The abundances 
of lithium inferred from observations of the very metal-poor halo stars (one- 
thousandth of solar and even lower) require almost no correction for chemical 
evolution. On the other hand, the status of helium-3 is in contrast to that of 
the other light elements. Although all prestellar D is converted to 3 He during 
pre-main sequence evolution, 3 He is burned to 4 He and beyond in the hotter 
interiors of most stars, while it survives in the cooler exteriors. For lower mass 
stars a greater fraction of the prestellar 3 He is expected to survive and, indeed, 
incomplete burning leads to the buildup of 3 He in the interior which may, or 
may not, survive to be returned to the interstellar medium [p3[ . In fact, some 
planetary nebulae have been observed to be highly enriched in 3 He, with abun- 
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dances 3 He/H ~ 10~ 3 ]24| |. Although such high abundances are expected in the 
remnants of low mass stars [^3[ ^5) , if all stars in the low mass range produced 
comparable abundances, we would expect solar and present ISM abundances of 
3 He to greatly exceed their observed values [^5[ |2(| ^?J. It is therefore necessary 
that at least some low mass stars are net destroyers of 3 He. For example, there 
could be "extra" mixing below the convection zone in these stars when they are 
on the red giant branch |2^, |3(| . Given such possible complicated histories of 
survival, destruction, and production, it is difficult to use the current Galactic 
and solar system data to infer (or, even bound) the primordial abundance of 
3 He. For this reason, we will not consider 3 He any further in this review. 

The generally favorable observational and evolutionary state of affairs for 
the nuclides produced during BBN is counterbalanced by the likely presence of 
systematic errors which are difficult to quantify and, in some cases, by a woefully 
limited data set. As a result, although cosmological abundance determinations 
have taken their place in the current "precision" era, it is far from clear that 
the present abundance determinations are "accurate" . Thus, the usual caveat 
emptor applies to any conclusions drawn from the comparison between the pre- 
dictions and the data. With this caution in mind we survey the current status 
of the data to infer "reasonable" ranges for the primordial abundances of the 
key light elements. 

3.1 Deuterium 

Deuterium is the ideal baryometer. As we have noted above the BBN-predicted 
D/H ratio is a strong function of the baryon-to-photon ratio r). A determination 
of the primordial abundance to 10%, leads to an 77 determination accurate to ~ 
6%. Furthermore, since deuterium is burned away whenever it is cycled through 
stars, and there are no astrophysical sites capable of producing deuterium in 
anywhere near its observed abundance ]3l| , any observed D-abundance provides 
a lower bound to its primordial abundance. Thus, without having to correct 
for Galactic evolution, the deuterium abundance inferred from UV observations 
of the local interstellar medium (LISM) {$2), D/H = (1.5 ± 0.1) x 10~ 5 (unless 
otherwise noted, observational errors are quoted at la), bounds the primordial 
abundance from below and the baryon-to-photon ratio from above p3fl . This 
value represents an average along 12 lines of sight in the LISM. Although they 
are not directly relevant to BBN, it is interesting to note that there have been 
several reports |34| [35| of a dispersion in ISM D/H abundances. It is not clear 
whether such variations are related to those inferred for the 3 He/H abundances 
in Galactic H 11 regions Q . 

Solar system observations of 3 He permit an indirect determination of the 
pre-solar system deuterium abundance (Geiss & Reeves 1972). This estimate 
of the Galactic abundance some 4.5 Gyr ago, D/H = (2.1 ± 0.5) x 10~ 5 (Geiss 
& Gloeckler 1998), while having larger uncertainty, is consistent with the LISM 
value. There has also been a recent measurement of deuterium in the atmosphere 
of Jupiter using the Galileo Probe Mass Spectrometer J3lj, which finds D/H 
= (2.6 ±0.7) x 10~ 5 . 
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To further exploit the solar system and/or LISM deuterium determinations 
to constrain/estimate the primordial abundance would require corrections for 
the Galactic evolution of D. Although the simplicity of the evolution of deu- 
terium (only destroyed) suggests that such correction might be very nearly in- 
dependent of the details of specific chemical evolution models, large differences 
remain between different estimates |3^, [59|. It is therefore fortunate that data 
on D/H in high-redshift, low-metallicity Lyman-a absorbers has become avail- 
able in recent years po|-p5|. It is expected that such systems still retain their 
original, primordial deuterium, undiluted by the deuterium-depleted debris of 
any significant stellar evolution. That's the good news. The bad news is that, at 
present, there are D-abundance determinations claimed for only four such sys- 
tems and that the abundances inferred for two of them appear to be inconsistent 
with the abundances determined in the other two. Here is a prime example of 
"precise" but possibly inaccurate cosmological data. There is a serious obstacle 
inherent to using absorption spectra to measure the deuterium abundance since 
the isotope-shifted deuterium lines are indistinguishable from velocity-shifted 
hydrogen. Such "interlopers" may have been responsible for some of the early 
claims [|o| of a "high" deuterium abundance |Q. Data reduction errors may 
have been the source of another putative high-D system. At present it seems 
that only three good candidates for nearly primordial deuterium have emerged 
from ground- and space-based observations. 

The absorption system at z = 3.572 towards Q1937-1009 was first studied 
by Tytler, Fan & Buries |lj who derived a low D/H = (2.3 ± 0.3 ± 0.3) x 
10~ 5 . Since an uncertain hydrogen column density, due to the saturated Lyman 
series profiles, was the largest source of uncertainty p7| , new, high quality, low- 
resolution spectra were obtained p2] ] , which, along with a new fitting procedure 
led to a revised abundance: D/H = (3.3 ± 0.3) x 10~ 5 ; notice the rather poor 
overlap with the original abundance. The z — 2.504 absorption system towards 
Q1009+2956 provides another potentially accurate D-abundance determination 
[|| D/H = (4.0 ± 0.7) x 10~ 5 . There are two other systems studied by Buries 
& Tytler (1998) whose derived D-abundances are consistent with these two, but 
whose uncertainties are much larger. The weighted mean of the two accurate D- 
abundance determinations leads to a 95% confidence range: 2.9 x 10~ 5 < D/H < 
4.0 x 10~ 5 . We adopt this range in our comparisons with the BBN predictions. 

We note that Levshakov, Kegel & Takahara ( jf8) , LKT) have used the data in 
[ pT| for the z = 3.572 system towards Q1937-1009, but with a different model for 
the velocity distribution of the absorbing gas, to derive a (95% confidence) range 
3.5 xl0~ 5 < D/H < 5.2 x 10~ 5 , which argues for a slightly higher abundance 
than suggested by the Buries & Tytler range. These same authors also 
used their model to reanalyze the Buries & Tytler |l3| data for Q1009+2956 
[|9) and they derive a 95% estimated range of 2.9 xl0~ 5 < D/H < 4.6 x 10~ 5 , 
now in excellent agreement with the Buries & Tytler value for this system. 
Recently, Levshakov, Tytler & Buries (]50|, LBT) have joined forces to apply 
this different model to a reanalysis of the z — 2.504 absorption system towards 
Q1009+2956, finding a consistent but slightly higher range (68%): D/H ~ (3.5- 
5.0) x 10~ 5 . 
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Although deuterium in the two high-redshift absorbers is consistent with a 
primordial abundance in the range 2.9 xl0~ 5 < D/H < 4.0 x 10~ 5 (or slightly 
higher accounting for the LKT and LBT analyses of the same data) , the deu- 
terium abundance derived for the one low-redshift absorber, the z = 0.701 
system towards Q1718+4807 observed with the GHRS on HST is significantly 
different. This data was first analyzed by Webb et al. Q who derived a 
very high deuterium abundance: D/H = (20 ± 5) x 10~ 5 . In contrast, LKT 
[ |5l| using the same data but their model for the velocity distribution of the 
absorbing gas, derive an abundance closer to those for the high-redshift ab- 
sorbers: 4.1 xl0~ 5 < D/H < 4.7 x 10~ 5 . Recently, Tytler et al. f|| use 
new Keck spectra to supplement the data from HST to derive a 95% range: 8 
xl0~ 5 < D/H < 57 x 10~ 5 , consistent with the Webb et al. (1997) estimate. 
Clearly the high-D abundance inferred from some analyses of this system are 
inconsistent with the low-D abundances derived from the other two, highcr- 
rcdshift systems. The sense of the discrepancy is puzzling since it is expected 
that the deuterium abundance should only decrease with time (decreasing red- 
shift). If, in fact, the high abundance is representative of the primordial value, 
then the other two absorbers should consist of gas most of which has been cy- 
cled through stars. The high redshifts and low metallicities of these systems 
suggest this is unlikely. If high D-abundances at high-z and low-Z are common, 
many systems like Q1718+4807 should present themselves for analysis. Tytler 
(1998) has argued that the absence (so far) of very many possible candidates 
suggests that either the abundance determination in the Q1718+4807 absorber 
is unreliable, or the Q1718+4807 absorber is anomalous. 

In anticipation of new data which may resolve this conundrum, we prefer to 
keep our options open and discuss the consequences of either of two (mutually 
exclusive) possibilities. For the low-D case we use the two high-z systems and 
adopt the Burles-Tytler 95% range: 2.9 xlO" 5 < D/H < 4.0 x 10~ 5 . For the 
high-D case we adopt the range: 1 xl0~ 4 < D/H < 3 x 10 based on the 
2cr range of Webb et al. [H. With account for the uncertainty in the BBN- 
predicted D-abundance at fixed 77, the lower bound to primordial D/H for the 
low-D case leads to an upper bound to 77 of: 7710 < 6.3, while the upper bound 
on D/H leads to a lower bound on 77 of: 7710 > 4.2. For the high-D option, the 
corresponding range in 77 is: 1.2 < 7710 < 2.8. In making these estimates we have 
been "overly generous" in the sense that the 77 values correspond to the "2a" 
uncertainties in the observational data and the "2cr" uncertainties in the BBN 
predictions. 

3.2 Helium-4 

As the second most abundant nuclide in the Universe (after hydrogen), the 
abundance of 4 He can be determined to high accuracy at sites throughout the 
Universe. To minimize the uncertainty inherent in any correction for the debris 
of stellar evolution, it is sensible to concentrate on the data from low-metallicity, 
extragalactic H11 regions |5^|-|59|]. Since each data set contains of order 40 re- 
gions, various analyses achieve statistical uncertainties in their estimate of the 
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primordial helium mass fraction < 0.003 (or, < 1%). Further, since the most 
metal-poor of these regions have metallicities of order 1/50 - 1/30 of solar, 
the extrapolation from the lowest metallicity regions to truly primordial intro- 
duces an uncertainty no larger than the statistical error. Although 4 He has 
already entered the era of "precision cosmology" , difficult to constrain system- 
atic uncertainties dominate the error budget. For example, using published 
data for 40 low- metallicity regions (excluding the suspect NW region of IZwl8), 
Olive & Steigman (OS) f| find: Y P = 0.234 ± 0.003 based on the data in 
| j52[ p3| . In contrast, from an independent data set of 45 low- metallicity re- 
gions with only slight overlap with that of OS, Izotov & Thuan (IT) |57) infer 
Yp = 0.244 ± 0.002. Clearly these results are statistically inconsistent. Several 
contributions to this discrepancy can be identified. Since the intensity of the 
helium recombination emission lines can be enhanced by collisional excitation 
[ |58| , corrections for collisional excitation are mandatory. In (54|, |5?| an attempt 
was made to use helium-line data alone (5 lines) to make this correction, in 
contrast to the traditional approach using information on the electron density 
derived from non-helium line data (see Skillman, Terlevich & Terlevich [ |60| for a 
discussion) . It is of great value that Izotov et al. |54| (ITL) and IT also analyze 
their data according to the traditional approach since this permits an estimate 
of the effect of their approach on the inferred primordial abundance. Using 
their data for 44 regions analyzed similarly to the data employed in OS, they 
would have derived Yp = 0.241 ± 0.002, reducing the discrepancy between OS 
and IT. Another source of systematic difference between the two analyses can 
be identified. By relying on helium (and hydrogen) recombination lines, any 
neutral helium (or hydrogen) present in the H n regions is invisible and must be 
corrected for. Since any such correction will be model dependent and uncertain, 
Pagel et al. |5^| restricted their attention to H n regions of "high excitation" for 
which this correction should be minimized. As a result they (and most of the 
data utilized by OS) make no ionization correction. In contrast ITL, through 
a misreading of published models of H n regions, make a correction for neutral 
helium while ignoring the (predicted) larger correction for neutral hydrogen in 
regions ionized by hot stars (metal-poor stars are hotter than the correspond- 
ing solar metallicity stars). Skillman, Terlevich & Terlevich |Q estimate the 
size of this correction to be of order 1% (AYp ks —0.002), further reducing the 
discrepancy between the IT and OS Yp estimates to « 0.005 rather than the 
original 0.010. Although IT eliminate the erroneous ionization correction from 
ITL in their more recent work, they actually derive a higher helium abundance. 
IT remark that this may be due to the higher temperatures in their new regions 
(compared to the ITL data set). 

At present potentially the most significant systematic uncertainty affecting 
the derived primordial abundance of helium appears to be that due to possible 
underlying stellar absorption (ITL; IT; Skillman, Terlevich & Terlevich (6^]). It 
has become clear that the helium abundance determination in the NW region 
of IZwl8 is likely contaminated by such absorption, resulting in an underesti- 
mate of the true abundance. Other regions in the OS and Olive, Skillman & 
Steigman fl56| (OSS) data sets may suffer similar contamination, biasing their 
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estimate of the primordial helium abundance to values which may be too low. 
In contrast, ITL/IT select their regions on the basis of the strength of the he- 
lium lines, avoiding those weak-lined regions which may be contaminated by 
underlying stellar absorption. If, indeed, they have been successful in avoiding 
this systematic error, their higher abundance estimate may be closer to the true 
value. But, through such selection they have run the risk of introducing a bias 
against finding low helium abundances. 

It is clearly crucial that high priority be assigned to using the H u region ob- 
servations themselves to estimate/avoid the systematic errors due to underlying 
stellar absorption, to collisional excitation and, to corrections for neutral helium 
and/or hydrogen. Until then, the error budget for Yp is likely dominated by 
systematic rather than statistical uncertainties and it is difficult to decide be- 
tween OS (and OSS) and IT. When account is taken of systematic uncertainties, 
they may in fact be consistent with each other. Therefore, in what follows, we 
will will adopt a generous "95%" range of 0.228 < Y P < 0.248 (cf, g). 

3.3 Lithium-7 

Cosmologically interesting lithium is observed in the Popn halo stars |fsi|-[|64f 
which are so metal-poor they provide a sample of more nearly primordial mate- 
rial than anything observed anywhere else in the Universe; the most metal-poor 
stars have less than one-thousandth the solar metallicity. Of course these halo 
stars are the oldest stars in the Galaxy and, as such, have had the most time to 
modify their surface abundances. So, although any correction for evolution mod- 
ifying the lithium abundance may be smaller than the statistical uncertainties 
of a given measurement, the systematic uncertainty associated with the dilution 
and/or destruction of surface lithium in these very old stars could dominate the 
error budget. Additional errors are associated with the modeling of the surface 
layers of these cool, low-metallicity, low-mass stars, such as those connected 
with stellar atmosphere models and the temperature scale. It is also possible 
that some of the observed Li is non-primordial, (e.g., that some of the observed 
Li may have been produced by spallation or fusion in cosmic-ray collisions with 
gas in the ISM f6§ 

There now exists a very large data set of lithium abundances measured in 
the warmer (T > 5800A"), metal-poor ([Fe/H] < —1.3) halo stars. Within the 
errors, these abundances define a plateau (the "Spite-plateau" ) in the lithium 
abundance - metallicity plane. Depending on the choice of stellar-temperature 
scale and model atmosphere the abundance level of the plateau is: A(Li) = 
12+log(Li/H)= 2.2±0.1, with very little intrinsic dispersion around this plateau 
value (e.g., [p2[). This small dispersion provides an important constraint on 
models which attempt to connect the present surface lithium abundances in 
these stars to the original lithium abundance in the gas out of which these 
stars were formed some 10 - 15 Gyr ago. "Standard" (i.e., non-rotating) stellar 
models predict almost no lithium depletion and, therefore, almost no dispersion 
about the Spite-plateau ]67|| . 

Early work on mixing in models of rotating stars was very uncertain, predict- 



12 



ing as much as an order of magnitude 7 Li depletion. Recently, Pinsonneault et 
al. |68| ] , building on progress in the study of the angular momentum evolution of 
low-mass stars |39| , constructed stellar models which reproduce the angular mo- 
mentum evolution observed for low-mass open cluster stars, and have applied 
these models, normalized to the open cluster data and to the observed solar 
lithium depletion, to the study of lithium depletion in main sequence halo stars. 
Using the distribution of initial angular momenta inferred from young open 
clusters for the halo stars leads to a well-defined lithium plateau with modest 
scatter and a small population of "outliers" (overdepleted stars) which is con- 
sistent with the data. Consistency with the solar lithium, with the open cluster 
stars, and with the (small) dispersion in the Spite-plateau may be achieved for 
depletion factors between 0.2 dex and 0.4 dex Q. 

The amount of depletion can also be limited |7(], [n], [72| by observations 
of 6 Li j73|. If the original 6 Li in halo stars is assumed to be as high as the 
solar value, an upper bound of 0.4 dex 7 Li depletion in rotational models is 
obtained from 6 Li data |}8|. Recent analysis suggests a more stringent 
(albeit model dependent) Li depletion limit of 0.2 dex based on constraints on 
the low metallicity ([Fe/H] w —2.3) production of 6 Li. Clearly 6 Li plays a vital 
role when it comes to constraining 7 Li depletion - the key issue to be resolved 
is the evolution of 6 Li in low metallicity environments and the data required are 
the simultaneous observations of the isotopes of Li, Be, and B in low metallicity 
halo stars. 

Very recently, Ryan, Norris & Beers (RNB) |64| have presented data for 23 
very metal-poor ([Fe/H] < — 2.5) field turnoff stars, chosen specifically to lie 
in a limited range of metallicity so as to facilitate the study of the dispersion 
in the Spite plateau. Although the limited data set subjects any conclusions 
to the uncertainties due to small number statistics, these data confirm previous 
suggestions |32j that there is very little dispersion about the plateau abundance. 
RNB claim evidence for a slope in the A(Li) vs [Fe/H] data (an increase of Li 
with Fe) . If real, this suggests that not all of the inferred lithium is primordial. 
In a recent analysis j74|, it is argued that 0.04 - 0.2 dex of the observed A (Li) 
could be post-primordial in origin. On the basis of the very small residual 
dispersion after accounting for the trend in A(Li) with [Fe/H], and with some 
"outliers" removed, RNB argue that their data (which may be statistics limited) 
is consistent with no dispersion and for an upper limit on the lithium depletion of 
0.1 dex. As discussed in the fraction of "outliers" is crucial for constraining 
rotationally mixed models. As of this writing most, if not all, evidence points 
to a rather limited depletion of no more than 0.2 dex, either in standard stellar 
models or in those including rotation. 

To err on the side of caution, we adopt a central value for the plateau 
abundance of A(Li) = 2.2 and we choose a ~ 2a range of ±0.1 dex so that our 
adopted "95%" range is: 2.1 < A(Li) < 2.3. If depletion is absent, this range 
is consistent with the lithium "valley" . For depletion > 0.2 dex, the consistent 
lithium abundances bifurcate and move up the "foothills" , although a non- 
negligible contribution from post-primordial lithium could move the primordial 
abundance back down again. 
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4 Confrontation Of BBN Predictions And 
Observational Data 

In the context of the "standard" model (three families of light or massless, 
two-component neutrinos), the predictions of BBN depend on only one free 
parameter, the nucleon-to-photon ratio ?/. Recalling that for To = 2.728 K, 
7/10 = 273Qb}i 2 , the baryon inventory of Persic & Salucci j75| may be used 
to set a very conservative lower bound, rjio > 0.25. From constraints on the 
total mass density and the Hubble parameter, the extreme upper bound on 
r) could be nearly three orders of magnitude larger. Over this large range in 
cosmologically "interesting" nucleon abundance, the predicted abundance of 
deuterium changes by more than eight orders of magnitude, from more than 
several parts in 10 3 to less than a part in 10 11 as can be seen in Figure 4, 
where the BBN predictions are shown over a wide range in rj. Over this same 
range in nucleon abundance, the lithium abundance varies from a minimum 
around 10 -10 to a maximum some two orders of magnitude larger, while the 
predicted primordial helium mass fraction is anchored between 0.2 and 0.3. 
Even the 3 He abundance, which we have set aside due to its uncertain Galactic 
evolution, varies from much higher than observed (> 10 -4 ) to much less than 
observed (sa 10~ 6 ). The key test of the standard, hot, big bang cosmology is 
to ask if there exists a unique value (range) of rj for which the predictions of 
the primordial abundances are consistent with the light element abundances 
inferred from the observational data. Since we have allowed for the possibility 
that one of the two current estimates of primordial deuterium from extragalactic, 
absorption studies could reflect the true abundance of primordial deuterium, our 
test must be done in two parts. Monte Carlo techniques have proven to be a 
useful tool in the analysis of the concordance between the BBN predictions and 
the observationally determined abundances of the light elements || Q |7^|-||79|. 
However, since for the purpose of this review we have taken a broad brush 
approach to the observational data, we limit ourselves to a simpler, more semi- 
quantitative discussion of this comparison. 

4.1 Low Deuterium 

From the two well-observed high-redshift absorption-line systems, we have adopted 
the Buries and Tytler 95% estimate for the primordial-D abundance: (D/H)p = 
2.9 — 4.0 x 10~ 5 . With allowance for the (2a) uncertainties in the BBN-predicted 
abundance (see Fig. 2), the consistent range of 77 is quite narrow: 7710 = 4.2 — 6.3. 
For this range in nucleon-to-photon ratio, the primordial lithium abundance is 
predicted (with account for the 2a uncertainties in the prediction) to lie in the 
range: A(Li)eBN = 2.1 — 2.8. In our discussion of the status of the lithium 
observational data we identified a range for its primordial abundance which has 
significant overlap with this predicted range (see Fig. 3): A(Li)p = 2.1 — 2.3. 
Thus, the D-constrained range of 7710 = 4.2 — 6.3, is consistent with the inferred 
primordial abundance of lithium, even allowing for ~ 0.2 dex stellar destruction 
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Figure 4: The predicted abundances as a function of 77. 



and/or galactic production. So far, so good. What of primordial helium? Over 
this limited range in 77, the predicted helium mass fraction varies but little. 
With account for the (small) uncertainty in the prediction (dominated for this 
range in 77 by the uncertainty in the neutron lifetime): Ybbn = 0.244 — 0.250. 
This range in the predicted primordial helium mass fraction, although on the 
high side, has significant overlap with the range inferred from observations of the 
low-metallicity, extragalactic H 11 regions: Yp = 0.228 — 0.248. For "low-D", the 
standard model passes this key cosmological test. For 7710 in the narrow range 
from 4.2 to 6.3, the predicted and observed abundances of deuterium, helium-4 
and lithium- 7 are in agreement (and, the predicted abundance of helium-3 is 
consistent with the abundances inferred for the interstellar medium and in the 
presolar nebula). 
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4.2 High Deuterium 

If, instead, the high abundance of deuterium derived from HST and Keck obser- 
vations of one relatively low-redshift absorption-line system is truly representa- 
tive of the primordial deuterium abundance, a different range for the nucleon-to- 
photon ratio is identified: 7710 = 1.2 — 2.8 (see Fig. 2). The predicted primordial 
abundance of lithium for this range is A(Li) — 1.9 - 2.7 revealing virtually per- 
fect agreement with the abundance derived from the very metal-poor halo stars 
in the Spite plateau. Over this same 77 range, the predicted helium mass frac- 
tion varies from Ybbn = 0.225 to 0.241. Here, too, the prediction is in excellent 
agreement with the observed abundance range. Thus, for "high-D" as well, the 
standard model passes this key cosmological test. 

4.3 Consistency With Non-BBN Estimates? 

Having established the internal consistency of primordial nucleosynthesis in the 
standard model, it is necessary to proceed to the next key test. Does the nucleon 
abundance inferred from processes which occurred during the first thousand 
seconds of the evolution of the Universe agree with estimates/bounds to the 
nucleon density in the present Universe? 

It is a daunting task to attempt to inventory the baryons in the Universe. 
Since many (most?) baryons may be "dark" , such approaches can best set 
lower bounds to the present ratio of baryons-to-photons. In their inventory of 
visible baryons, Persic & Salucci |7q ] estimate for the baryon density parameter: 
f^B ~ 0.0022 + 0.0015/15Q 1 ' 3 , where /150 is the Hubble parameter in units of 
50 km/scc/Mpc. For a lower bound of Ho > 50 km/sec/Mpc, this corresponds to 
a lower bound on 77 of: 7710 > 0.25, entirely consistent with our BBN estimates. 
More recently, Fukugita, Hogan & Peebles [^| have revisited this question. 
With subjective, but conservative estimates of the uncertainties, their lower 
bound to the global budget of baryons (for Ho > 50 km/sec/Mpc) corresponds 
to a much higher lower bound: 7710 > 1.5, which is still consistent with the 
"low- 77" range we identified using the high D results. A possible challenge to 
the "low-77" case comes from the analysis of Steigman, Hata & Felten |[l| who 
used observational constraints on the Hubble parameter, the age of the Universe, 
the "shape" parameter, and the X-ray cluster gas fraction to provide non-BBN 
constraints on the present density of baryons, finding that 7710 > 5 may be 
favored over 7710 < 2. Even so, a significant low-77, high-D range still survives. 

5 Constraints from BBN 

Limits on physics beyond the standard model are mostly sensitive to the bounds 
imposed on the 4 He abundance. As described earlier, the 4 He abundance is pre- 
dominantly determined by the neutron-to-proton ratio just prior to nucleosyn- 
thesis; this latter is set by the competition between the weak interaction rates 
and the universal expansion rate. Modulo the occasional free neutron decay, the 
neutron-to-proton ratio "freezes-out" at a temperature ~ 800 keV. While the 
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weak interaction rates converting neutrons to protons and vice- versa are "fixed" , 
there may be room for uncertainty in the expansion rate which depends on the 
total mass-energy density. For example, the presence of additional neutrino fla- 
vors (or of any other particles which would contribute significantly to the total 
energy density) at the time of nucleosynthesis would increase the total energy 
density of the Universe, thus increasing the expansion rate, leading to an earlier 
freeze-out, when the temperature and the n/p ratio are higher. With more neu- 
trons available, more 4 He can be synthesized. In the standard model the energy 
density at a temperature of order 1 MeV is dominated by the contributions from 
photons, electron-positron pairs and three flavors of light neutrinos. We may 
compare the total energy density that in photons alone through N which counts 
the equivalent number of relativistic degrees of freedom. 

p = (N/2) Pl (4) 

In the standard model at T ~ 1 MeV, Nsm = 43/4, so that we may account for 
additional degrees of freedom by comparing their contribution to p to that of 
an additional light neutrino species 

N = N SM + 7/8AN, (5) 

For AN U sufficiently small, the predicted primordial helium abundance scales 
nearly linearly with AN V : AY ks 0.013A7V,,. Hence, any constraints on Y 
lead directly to bounds on AN U [ p2[ . However, it is worth recalling that the 
constraint is, ultimately, on the ratio of the Hubble parameter (expansion rate) 
and the weak interaction rate at BBN, so that changes in the weak and/or 
gravitational coupling constants can be similarly constrained Jll], |8j| . Here we 
will restrict our attention to the limits on N v and on neutrino masses from BBN. 
Although likelihood methods have been used to obtain more exact limits on N v 
[p3, again here we adopt a simpler, more broad brush approach. Many of the 
limits on particle properties were recently reviewed in Q ■ 

Given the observational upper bound on Yp of 0.248 and a predicted lower 
bound of 0.244 (for low-D), there is room for an increase in the BBN-predicted 
4 He of AY = 0.004. From the scaling of Y with AN V , we derive an upper limit 
to AN V of AN V < 0.3. It should be cautioned that this bound is really less 
stringent than a true "2ct" upper limit, since we have chosen 2a ranges both 
in the predicted and the observed deuterium and helium abundances. Even so, 
for low-D this constraint is already good enough to permit an exclusion of any 
"new", light scalars (which would count as AN V = 0.57), as well as a fourth 
neutrino. For high-D we predict a lower bound of Y = 0.225 to be compared 
with the observed upper bound of with Y = 0.248, and using the same argument, 
we derive an upper bound of AN V < 1.8. 

It should be noted that the limit derived above is not restricted to full 
strength weak interaction neutrinos. In fact, since we know that there are 
only three standard neutrinos, the limit is most usefully applied to additional 
particle degrees of freedom which do not couple to the Z°. For very weakly 
interacting particles which decouple very early, the reduced abundance of these 
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particles at the time of nucleosynthesis must be taken into account pq| . For a 
new particle, x, which decoupled at Td > 1 MeV, conservation of entropy relate 
the temperature of the xs to the photon/neutrino temperature (T) at 1 MeV, 
((T x /T) 3 = ((43/47V(T rf ))). Given g B{F) boson (fermion) degrees of freedom, 

^^Ef^ + Ef^ 4 - (6) 

As an example of the strength of this bound, models with right-handed interac- 
tions, and three right-handed neutrinos, can be severely constrained since the 
right-handed states must have decoupled early enough to ensure that 3(T l/R /T„ L ) 4 
< AN U . Using the high D limit to N„, three right-handed neutrinos requires 
N(Td) > 15, implying that Td > 40 MeV. In contrast, the low D limit requires 
that N(Td) > 60 so that Td > 300 MeV. If right-handed neutrino interactions are 
mediated by additional gauge interactions, associated with some scale Mz', and 
if the right handed cross sections scale as M^/ 4 , then the decoupling temperature 
of the right handed interactions is related to Mz 1 by (T^/X^) 3 oc {Mz 1 /Mz) 4 
which, for TdL ~ 3 MeV requires TdR ^ 40(300) MeV, the associated mass scale 
becomes Mz' > 0.6(2.8) TeV! Note that this constraint is very sensitive to the 
BBN limit on N v . 

Many other constraints on particle properties can be related to the limit 
on N v . For example, neutrinos with MeV masses would also change the early 
expansion rate, and the effect of such a neutrino can be related to that of 
an equivalent number of light neutrinos [g7| ^8|, |9(|. A toy model which 
nicely contains ways to both increase and decrease 4 He production relative to 
standard BBN is the case of a massive v T [ p8[ . The two relevant parameters 
are the v r mass and lifetime. A v T which is stable on BBN timcscales (i.e., , 
T i> ft 100 sec) and has a mass greater than a few MeV will increase Yp relative 
to standard BBN. This is because such a neutrino still decouples when it is 
semi-relativistic and so its number density is comparable to that of a massless 
neutrino. However, its energy density at the onset of BBN is much greater 
than that of a massless neutrino since its mass is significantly greater than 
the temperature. Therefore, weak interactions decouple earlier, increasing the 
neutron-to-proton ratio at freeze out and thus the amount of 4 He. For example, 
a limit of N u < 4 translates into a mass limit on a relatively stable (on the 
time-scale of BBN) neutrino of m v < 0.4 MeV for a Dirac-mass neutrino and 
m„ < 0.9 MeV for a Majorana-mass neutrino J9(J. Just the opposite can occur 
if such a v r decays rapidly compared to BBN timcscales. The rapid decays and 
inverse decays keep the v T s in equilibrium much longer than do the conventional 
weak interactions so that their number density, along with their energy density, 
is exponentially suppressed. A typical example is a relative decrease in Yp of 
about 0.01 for a v T with a mass of ~ 10 MeV and a lifetime [y T — > -\- cf) where 
is a Majoron) of 0.1 sec. 
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6 Conclusions 



In this precision era of Cosmology the BBN abundances are predicted with great 
accuracy in the standard model. The statistical uncertainties in the primordial 
abundances of the light nuclides inferred from the observational data are also 
very small. However, there is evidence that the derived abundances may be 
subject to systematic errors much larger than the statistical errors. This is 
particularly evident for deuterium where the D/H ratio derived for two, low 
metallicity, high redshift absorption systems differs by a factor of 5 - 10 from 
that inferred for a third such system. For 4 He, two determinations of the pri- 
mordial mass fraction differ from each other by 2 - 3 times the statistical error. 
Their differences may be traced to differing treatments of the corrections for 
collisional excitation and ionization and the data sets may be contaminated by 
some cases of underlying stellar absorption. Although a clear, accurately deter- 
mined "plateau" is evident in the Li vs. Fe relation for the metal-poor halo stars, 
the level of the plateau is subject to uncertainties in the metal-poor star tem- 
perature scale and atmosphere models. In addition, there may be non-negligible 
corrections (larger than the statistical uncertainties) due to depletion of surface 
lithium in these very old stars, as well as enhancement due to post-BBN produc- 
tion. Nonetheless, despite these nagging uncertainties, the agreement between 
the predictions of standard BBN and the observed abundances is impressive. 
The standard model passes this key test with flying colors. 

Given the dichotomy in the possible primordial abundance of deuterium, we 
have considered two possibilities. For the "low-D" option, we identify a "high- 
rf range (at 95% confidence): rj W w 4.2 — 6.3. In this range the predicted 
abundances of 3 He, 4 He and 7 Li are consistent with the primordial abundances 
inferred from observations (see Figures 1-3). For rj in this range the baryon 
density parameter is restricted to: Q,^h? w 0.015 — 0.023 which, for H = 
70 km/s/Mpc corresponds to: f^B ~ 0.03 — 0.05. Using the upper bound to 
Yp from the data along with the lower bound to ij leads to a "high-77" bound to 
the number of "equivalent" light neutrinos: < 3.3. For the "high-D" option 
a "low-77 range is identified: ijiq w 1.2 — 2.8. In this range as well there is overlap 
between the predicted and observed primordial abundances of 3 He, 4 He and 7 Li. 
For this "low-77" range, n B h 2 w 0.004 - 0.010 which, for H = 70 km/s/Mpc 
corresponds to: SIb ~ 0.01 — 0.02. In this range the upper bound to the number 
of equivalent light neutrinos is much less restrictive: N„ < 4.8. As a key probe 
of early Universe Cosmology and of particle physics (standard model as well as 
beyond the standard model), BBN is alive and well. 
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